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Abstract. BPM 37093 is the only hydrogen-atmosphere white dwarf currently known which has sufficient mass (~1.1 M Q ) to 
theoretically crystallize while still inside the ZZ Ceti instability strip (T e g ~ 12,000 K). As a consequence, this star represents 
our first opportunity to test crystallization theory directly. If the core is substantially crystallized, then the inner boundary for 
each pulsation mode will be located at the top of the solid core rather than at the center of the star, affecting mainly the average 
period spacing. This is distinct from the "mode trapping" caused by the stratified surface layers, which modifies the pulsation 
periods more selectively. In this paper we report on Whole Earth Telescope observations of BPM 37093 obtained in 1998 and 
1999. Based on a simple analysis of the average period spacing we conclude that a large fraction of the total stellar mass is 
likely to be crystallized. 
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Kanaan et al.: WET Observations of BPM 37093 



1. Introduction 



Table 1. Journal of observations in 1998 (XCOV 16). 



Since 1960 most astronomers have agr eed that cool w hite 
dwarfs must eventual ly crystallize ( Kirzhn it J 1 96fjt lAbrikosovl 
1960; Salpeter 1961). The process theoretically begins when 
the electrostatic interaction between the ions becomes much 
larger than the thermal energy. This effect is based on such well 
known physics that it has become widely accepted without ever 
having been tested empirically. 



Run 


Telescope 


Date 


Start 


Length 








(UT) 


(s) 



BPM 37093 is a ZZ Ceti sta r dKanaan et all 19921) with an 
unusually high mass (1.10 M Q , Bergeron et alJl2004) . White 
dwarfs this massive are subject to much higher pressures and 
densities in their cores, and we expect a 1 .0 M white dwarf to 
begin crystallizi ng at temperatures within o r above the ZZ Ceti 
instability strip fWoo dill 9921: IWinget et allll 997h . Our goal in 
observing BPM 37 093 with the Whole Earth Telescope (WET, 
iNafheret all ll 990) was to obtain seismological data to deter- 
mine whether or not the stellar interior is crystallized. The fun- 
damental objective was simply to detect as many independent 
pulsation modes as possible, and then to compare the observed 
frequencies with those calculated from white dwarf models that 
have been artificially crystallized to various degrees. 

For years we have faced a troubling ambiguity between 
the effects of the crystallized mass fraction (M cr ), the hy- 
drogen layer mass (Mh), and the stellar mass (M») and 
effective temperature (T e ff). Changes to these four charac- 
teristics of white dwarf models can all modify the aver- 
age spacing between the c alculated pulsation periods (see 
[Montgomery & Wingetl Il999t eq. 7). Fortunately, the lat- 
ter two quantities can be constrained by fitting model at- 
mospheres to spectroscopic observations, but the others can 
only be determined through asteroseismology. Recent improve- 
ments in our theoretical description of the composition transi- 
tion zones between the stratified surface layers in our models 
dCorsico et alJ2002llAlthaus et alJ20"03l) have helped to reduce 
the degeneracy between Mh and M cr . However, the huge num- 
ber of possible parameter combinations, and the need for an 
efficient method of exploring them, remained serious obsta- 
cles to progress until recently (Metcalfe & Charbonneau 2003; 
Met calfe. Montgomery. & K anaan 20oJ 



Like the cool DAV G 29-38 jKleinman et alJ 1998), 
BPM 37093 exhibits irregular modulations in the amplitudes 
of its pulsation modes. On one occasion all of t he modes van- 
ished below the detection threshold of ~ 1 mma dKanaan et alJ 
119981) . However, the modes that disappeared were observed to 
reappear later with the same pulsation frequencies. This gives 
us confidence that we can learn more about this star by using 
the full set of frequencies that have been observed over time — 
a concept known as "time-ense mble" asteroseismology, pio- 
neered hv lKleinmane7aTl iri998). In this paper we report WET 
observations obtained in 1998 and 1999, and we use the iden- 
tified pulsation periods to define a range for the average period 
spacing. This allows us to constrain M» and M„ by following 
the analysis of Montgomery & Winget ( 1999) with an updated 
prescription for the envelope composition transition zones. 
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2. Observations 

In 1996 and 1997, observations of BPM 37093 were obtained 
from the 0.9 m telescope at CTIO and the 1.6 m telescope at 
Observatorio Pico dos Dias (OPD, Brazil) respectively. The 
two initial goals of these observations were: 1) to identify as 
many pulsation modes as possible, to help constrain asteroseis- 
mological model fitting, and 2) to find stable pulsation modes 
suitable for measuring the rate of period change (P), as has been 
done for other white dwarfs dCosta et all 1 19991: iKeoler et all 
2000al lMukadam et all2003l) . Further analysis (see sect.|3j re- 
vealed that no modes were stable enough to use for P mea- 
surements. After these two attempts to obtain single-site data 
on BPM 37093, it became clear that we would be unable to 
resolve the pulsation spectrum of this star from a single obser- 
vatory. By 1997 we had already accumulated more than 100 
hours of photometry on BPM 37093, which led to the identi- 
ficatio n of only 4 pulsatio n modes with highly variable ampli- 
tudes (Kanaan et al. 1998J). 



Kanaan et al.: WET Observations of BPM 37093 



Table 2. Journal of observations in 1999 (XCOV 17). 
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Fig. 1. Nightly FTs of BPM 37093 during ten nights of obser- 
vations from CTIO in 1996. The amplitudes are highly variable 
on short timescales. All of the modes vanish below the detec- 
tion threshold on the final night. 



provide an independent determination of the spherical degree 
{£) for each pulsation mode. Unfortunately, this method could 
be applied only to the modes with the highest amplitudes, and 
the observations were not sensitive enough to distinguish be- 
tween £ — 1 and t — 2. In this paper we infer £ based only on 
the average period spacing between modes. 



BPM 37093 was chosen as the southern primary target for 
a Whole Earth Telescope campaign (XCOV 16) in 1998, and 
again in 1999 (XCOV 17) to coincide with simultaneous obser- 
vations using the Hubble Space Telescope (HST). A journal of 
observations for the data obtained for these two campaigns is 
shown in Tables ^andE] Overall, we obtained more than 142 
hours of data in April-May 1998 with a duty cycle of 50% dur- 
ing the central 10 days, and an impressive 180 hours in April 
1999 with a better duty cycle of 65% during the central 10 
days. The latter observations included almost complete cover- 
age during the two scheduled HST visits near the middle of the 
campaign, and preliminary results were reported bv lNitta et all 
(2000)." 

The primary goal of the HST o bservations was t o use the 
limb darkening method devised by iRobinson et all d 19951) to 



3. Frequency Analysis 

One of the initial goals of our observations was to identify 
pulsation modes stable enough to measure the rate of period 
change (P). As a white dwarf star cools over time, the slow 
change in the thermal structure should lead to a detectable in- 
crease in the pulsation periods. We expect that the periods in 
a crystallizing star should change more slowly than in other 
ZZ Ceti stars because the associated release of latent heat 
causes the temperature to drop more slowly than in a non- 
crystallizing star. 

However, the observations from CTIO in 1996 made it clear 
that P measurements would not be possible for BPM 37093. In 
Fig.[2we show the evolution of daily Fourier Transforms (FTs), 
night by night during the 10 night run. On the tenth night, the 
pulsations vanished below the detection threshold. Comparing 
the other panels of Fig. ^ it is also clear that the amplitudes 
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Fig. 2. Fourier Transforms and window functions at the same scale for the Whole Earth Telescope observations of the ZZ Ceti 
star BPM 37093 obtained during a) the XCOV 16 campaign in 1998, and b) the XCOV 17 campaign in 1999. 



of the detected modes are highly variable. iKanaan etalJ ll998) 
demonstrated that these amplitude changes must be intrinsic to 
the star, rather than due to the beating of closely spaced modes, 
since there was no correlation between the changes in ampli- 
tude and phase. Such intrinsic behavior has been reporte d for 
other ZZ Ceti stars on a timescale of months jKleinman et alJ 
Il998l) . but not from night to night as observed in BPM 37093. 
These short-timescale variations will lead us to derive average 
amplitudes from the FTs of long data sets, but will still allow 
us to identify the pulsation modes which have the highest am- 
plitudes and the longest lifetimes. 

The FTs of the two WET campaigns on BPM 37093 are 
shown in Fig. |2 In each panel we have marked all of the 
peaks that exceed 4 times the observational noise level. During 
XCOV 16 we detected many more modes than were evident 
in the 1996 observations — of the 8 frequencies detected, only 
2 coincide with those seen previously. Comparing this with 
XCOV 17, it seems clear that the detection of more independent 
modes in XCOV 16 can be attributed to the star itself, rather 
than to differences in our detection threshold. During the sec- 
ond WET campaign on BPM 37093 we detected 7 frequencies, 
of which 4 coincide with others seen previously. Only the 512 s 
and 661 s modes are unique to the data obtained for XCOV 17. 
The full list of frequencies detected in these two WET cam- 
paig ns is shown in TableEl a long with the m ode identifications 
of M etcalfe. Montgomery. & Kanaanli 2004) 34 for those m odes 
found in the preliminary analysis of Kanaan et al. (2000). 

We can understand the erratic behavior of BPM 37093 by 
cons idering pre vious o bservations of other ZZ Ceti pulsators. 
iKleinman et alJ {l998) documented very similar results for the 
star G 29-38, and proposed that each set of observations can 
provide a subset of the full spectrum of normal modes. The 
fundamental idea behind this hypothesis is that the underly- 
ing structure of the star does not change on the timescales of 
the observed amplitude variations. Instead, something in the 



Table 3. Frequencies detected from two WET campaigns. 
XCOV Frequency Period (Amplitude) Identification 
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For the 582 s mode, we adopt the £ = 2 identification 



excitation mechanism must select certain modes and exclude 
others in a manner that varies over time. Setting aside the ques- 
tion of what specifically causes these amplitude modulations, 
we only need to assume that the observed frequencies are, in 
each case, normal modes of oscillation that can be described 
by spherical harm onics. Such an assumption rests on a firm 
body of evidence jRobinsonetaljri982t iKepler et alJ EoOOb: 
IClemens et alJl2000l) . and is also supported by the fact that 
modes which disappear below the detection threshold are ob- 
served to reappear later with the same frequencies . 

4. Interpretation 

We adopt for our analysis the full set of frequencies that 
have been observed in BPM 37093 over time, giving pref- 
erence to the WET campaigns for their superior frequency 
precision. For the isolated frequencies we assume that each 
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mode has an azimuthal order m = (see Metcalfe! [2003). 
and for the modes consisting of two closely-spaced frequen- 
cies we use the average of the two. For modes that were 
observed in both campaigns, we use the frequency from the 
observation with the highest amplitude. This leads to a total 
of 7 independent modes that have been id entified as t -2 
bv lMetcalfe. M ontgom ery. & Kanaanl J2004h . The average pe- 
riod spacing between consecutive radial overtones for these 7 
modes is (AP) = 17.6 + 1.1 s. This implies that the ( — 2 
identification should be correct for the majority of the modes, 
since an I — 1 identific ation would yield a much lar ger period 
spa cing ((AP) ~ 30 s.|Mo ntgomerv & Wing etll999l) . 

Montgomery & Wingetl Jl99 9) performed a detailed study 
of the various model parameters that could affect the average 
period spacing for 1-2 modes. They defined a scaling rela- 
tion for (AP) which had contributions from variations to four 
parameters, including the stellar mass and effective temperature 
(M*, r e ff) and the hydrogen layer mass and crystallized mass 
fraction (Mh, M a ). Fortunately, M* and P e ff can both be con- 
strained by spectral profile fittin g. The most rece nt estimates 
for BPM 37093 come from Bergeron et who found 

Af, = 1.10 M and T eff = 11 730 K. The values of the other 
paramete rs can only be dete rmined through asteroseismology, 
but lMontgomerv & Wingetl i 19991) found a troubling degener- 
acy between Mh and M cr that could not be broken using only 
the observed (AP). 

To get beyond the degeneracy, we need to use the individ- 
ual pulsation periods in addition to (AP). Fundamentally, this 
is possible because variations to Mh change the individual pe- 
riods through "mode trapping", while variations to M cr affect 
mainly the average period spacing. If BPM 37093 is partially 
crystallized, the inner boundary for each pulsation mode will 
be located at the top of the solid core rather than at the center 
of the star. This reduces the size of the resonant cavity, increas- 
ing the average period spacing and modifying the periods of all 
of the modes. By contrast, the hydrogen layer produces a sharp 
chemical gradient somewhere near the surface. Any individual 
mode whose eigenfunction is large in this region can interact 
with the chemical gradient and become "trapped" — its period 
will be shifted more than the periods of other modes. Clearly, 
each mode will be modified in a different manner through this 
int eraction. 

Montgomery & Wingetl \ 19991) understood that the signa- 
tures of Mh and M a are in principle distinct. However, by fo- 
cusing only on the average period spacing they were unable 
to place strong constraints on eit her Mh or M„. Our pr esent 
models are somewhat improved: Montgomery & Win get used 
hydrogen profiles that were derived assuming diffusive equilib- 
rium in the trace element approximation. This produced unre- 
alistically sharp chemical gradients at the base of the hydro- 
gen layer, leading to stro nger mo de trapping in their models. 
This was demonstrated bv lCorsico et alJll2002h . who compared 
models that assumed diffusive equilibrium in the trace element 
approximation with models that computed the abundance pro- 
files based on time-dependent diffusion calculati ons. In a recent 
extens ion of this work to massive ZZ Ceti stars. lAlthaus et all 
(2003) described an improved method of calculating diffusive 
equilibrium profiles that compare favorably with the fully time- 
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Fig. 3. The average period spacing of a small grid of models 
with various combinations of Mh and M cl . The lcr range of the 
observed average period spacing for BPM 37093 is shown as a 
shaded area, and the circled point indicates the model with the 
smallest rms difference between the observed and calculated 
periods (see text for details). 

dependent results (see their Fig. 18). We have incorporated this 
method of computing the hydrogen abundance profiles into the 
code used bv Montgo mery & Wingetl Jl999l) . However, since 
the sharpness of the hydrogen transition zone should mainly 
affect the mode trapping properties of the models, we expect 
that our new average pe riod spacings will differ only s lightly 
from those computed by Montgo mery & Wingetl Jl999l) . 

As a simple illustration of the potential of our observa- 
tions, we calculated (AP) for a small grid with various com- 
binations of Mh and M a . We fixed the mass, temperature, 
and helium layer thickness to the values used for Fig. 10b of 
Montgomery & Wing etl if 1999), but we assumed a uniform O 
core. We show this grid of models in Fig. [3] with the shaded l<x 
range of the average period spacing from the WET observa- 
tions of BPM 37093. As expected, the average period spacing 
of the 0% crystallized model is vir tually identical to that found 
bv Montgomerv & Winget ( 1999). However, due to the differ- 
ent assumed C/O profiles, the crystall ized curves have sh ifted 
with respect to the results of Montgomery & Wingetl lll999l) . 

Unfortunately, the degeneracy between Mh and M CI is still 
present, but we have not yet used the hidden third dimension of 
the grid: at each point we have also computed the root-mean- 
square differences between the observed and calculated peri- 
ods (cr P ), to exploit the information contained in the individual 
modes. Although the observed range of (AP) can accommodate 
values of log(M H /M„) between -5 (M a - = 80-90%) and -8 
(uncrystallized), not all hydrogen layer masses are equally suc- 
cessful at reproducing the individual modes. The model with 
log(Mtf/M„) = -6 and M cr = 50% has cr P = 1.08 s, which 
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is substantially better than anything else in this small grid (the 
next best model has <x P = 1 .70 s). A theoretical model with this 
same set of structural parameters is expected to be between 66- 
92% crystallized for a C/O mixture, or even mo re crystal lized 
if the core is composed of an O/Ne mixture (see lCorsico et alJ 
120041 for some recent calculations). 

Of course, the individual modes also contain information 
about the mass, temperature, helium layer, and core composi- 
tion. If we have fixed these parameters incorrectly, it is likely 
that we have found a locally optimal match to the observed 
periods rather than the global solution. Montgomerv & Winget 
(1999) recognized this difficulty, and discussed the need for an 
automated procedure to search this enormous parameter-space. 
Such a procedure, based on a parallel genetic algorithm, has 
recently been developed ( Metcalfe & C harbon neaull2003l) and 
applied to this problem. The initial results from this large-scale 
exploration of the models have been pub lished separately by 
iMetcalfe. Montgomery. & KanaanlJ20 04). who present a more 
detailed discussion of the initial model fitting results. 

5. Discussion 

The Whole Earth Telescope has once again deciphered the 
complex pulsation spectrum of an astrophysically interesting 
white dwarf star. BPM 37093 is the only known ZZ Ceti star 
massive enough to allow a seismological test of crystallization 
theory, and previous attempts to understand it from single site 
observations were not successful. The superior frequency pre- 
cision and extended coverage of two WET campaigns finally 
allowed us to document a series of 9 independent pulsation 
modes in this star. While it is always useful to search for ad- 
ditional frequencies to help constrain the model fitting, the ob- 
servational requirements of this project have now largely been 
satisfied. 

The limiting factor in our ability to test crystallization the- 
ory through asteroseismology of BPM 37093 is now computa- 
tional. The initial study of Metcalfe. Montgomerv. & Kanaan 
(2004) was limited to several fixed masses and core compo- 
sitions, but all of their optimal models suggested that a large 
fraction of the core is crystallized — a result that is qualita- 
tively supported by our simple analysis of the average pe- 
riod spacing. We have only fit M cr to the nearest 0.1 M t , but 
Montgomery & Wingetl i 19991) showed that the pulsations are 
sensitive to changes of 0.01 M, in this parameter. As comput- 
ers get faster, we should be able to extend the genetic algorithm 
based model fitting method to treat M a , M, and the composi- 
tion of the liquid portion of the core as fully adjustable param- 
eters. In the meantime, we should apply the same fitting pro- 
cedure to lower mass ZZ Ceti stars that are not expected to be 
crystallized, as a check of the method. 

The Sloan Digital Sky Survey (SDSS) has recently dis- 
covered several new ZZ Ceti stars th at may also be massiv e 
enough to test crystallization theory (Mukada m et al J 12004). 
All of the new pulsators in the SDSS sample are significantly 
fainter than the previously known white dwarfs, so larger tele- 
scopes will probably be required to resolve their pulsation 
spectra. Several cu rrent an d planned space m issions (MOST, 
IWalker et al.l2003t COROT. Ka olin et alJ2002T) promise to rev- 



olutionize the field of asteroseismology in the coming decade. 
However, none of them is likely to contribute to our under- 
standing of pulsating white dwarfs stars, which are too faint 
to be observed by these satellites. The Whole Earth Telescope 
remains the instrument of choice for these fainter targets. The 
availability of larger ground-based telescopes may allow the 
WET concept to be extended to 2 m and 4 m class instruments. 
With a growing harvest of new objects from SDSS, the future 
looks bright for white dwarf asteroseismology. 
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